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ABSTRACT 

A completely analytical theory of chirped-pulse oscillators is presented. The theory is based on an approximate 
integration of the generalized nonlinear complex Ginzburg-Landau equation. The obtained parametric space of 
a chirped-pulse oscillator allows easy tracing the characteristics of both solid-state and fiber oscillators operating 
in the positive dispersion regime. 
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1. INTRODUCTION 

In the last decade, femtosecond pulse technology has evolved rapidly and allowed achieving a few-optical-cycle 
pulse generation directly from an oscillator.^ Applications of such pulses range from medicine and micro- 
machining to fundamental physics of light-matter interaction at unprecedented intensity level and time scale. 
High-energy laser oscillators nowadays have reached the intensity level of the order of 10'^'* W/cm^, which 
allows high-intensity experiments such as direct gas ionization.® To achieve these regimes, about- and over- 
microjoule pulse energies are required. Such energy frontiers have become achievable due to the chirped pulse 
amplification outside an oscillator.^ However, the amplifier technology is i) complex, ii) expensive, iii) accessible 
pulse repetition rates lie within the kHz range, and iv) noise amplification is unavoidable. 

It is desirable to find a road to the direct over-microjoule femtosecond pulse generation at the MHz pulse 
repetition rates without an external amplification. To date, a promising approach has been proposed. It is based 
on a considerable decrease of the oscillator repetition rate. ^ The catch is that a long-cavity oscillator suffers from 
strong instabilities caused by enhanced nonlinear effects owing to increase of the pulse peak power P (0). The 
leverage is to stretch a pulse and, thereby, to decrease its peak power below the instability threshold. Recently, a 
critical milestones, demonstrating the feasibility of this approach, has been achieved for the Ti:sapphire oscillators 
operating both in the negative- (NDR)^' and positive-dispersion regimes (PDR),^^' the near-infrared Yb:YAG 
thin-disk oscillators operating in the NDR,^^"^^ and the fiber oscillators operating in the all-normal dispersion 
(ANDi) regime (that is the PDR by definition). 

The fundamental difference between the NDR and the PDR is that, in the first one, the Schrodinger soliton 
develops. Since the soliton peak power P (0) has to be lower than some threshold value Pth in order to avoid the 
soliton destabilization, the maximum reachable energy can be estimated as £' = 2PthT (T is the soliton width). 
That is the energy scaling requires the pulse stretching. However, the latter results from substantial growth of 
the group-delay dispersion (GDD) (quadratically with energy^*). As a result, the energy scaling requires a huge 
negative GDD, the obtained chirp-free soliton has a large width, and it is not compressible linearly. 

In the PDR,^^'^'^ the pulse is stretched and its peak power is reduced due to large chirp. ^^'^^ The chirp 
compensates the spectrum narrowing with energy and the pulse becomes to be compressible linearly down to 
T K, 2/A, where A is the spectrum half-width. The issue is that the chirped solitary pulse (CSP) is a dissipative 
soliton, that is it develops in a dissipative nonlinear system and, as a result, there is no a uniform description of 
its properties and dynamics, because the underlying nonlinear equation (so-called, nonlinear complex Ginzburg- 
Landau equation, CGLE) is not integrable.^'^^^"' 

In this work, I propose the approximate method of integration of the generalized nonlinear CGLE and 
demonstrate that the CSP is its solitary pulse solution with reduced dimension (2 or 3) of the parametric space. 
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As a result, the CSP characteristics become easily traceable on the two-dimensional diagram ("master diagram"). 
Comparison of the PDR parameters demonstrates that the CSPs formed in the ANDi fiber oscillator and in the 
CPO: i) lie within the distinct sectors of the unified master diagram, ii) belong to the distinct branches of solution, 
and iii) vary with parameters in different ways. Comparison of the models based on the different versions of 
the master equation is carried out. The phenomenon of concave spectrum is attributed to the quintic self-phase 
modulation. 



The nonlinear CGLE is the generalized form of the master mode-locking equation^^'^'^'^^ and provides an ade- 
quate description of mode- locked oscillators (both fiber and solid-state) . Its soliton-like solutions (or dissipative 
solitons) pattern the laser pulses. Such an approach is well-grounded if i) T ^ 2-k jtot^ (loq is the carrier frequency 
of laser field) and ii) relative variation of laser field during one cavity round-trip is small. 

Let u (2:, t) be a slowly-varying field amplitude, z be a propagation coordinate normalized to the cavity period, 
i be a local time. The generalized CGLE is 



where |m| is the instant power, a is the saturated net-loss, and a is the squared inverse spectral bandwidth of 
oscillator (as a rule, it is defined by gain bandwidth). Parameter /? is the net-GDD coefficient; 7 is the self-phase 
modulation (SPM) coefficient, and x describes a high-order correction to it (i.e. the quintic SPM). Function 

/ (l'"!^) models the self-amplitude modulation (SAM) in an oscillator and its form depends on the mode-locking 
mechanism. 

The CSP develops in the PDR under combined action of two mechanisms: the pure phase and dissipative 
ones. The first one results from a balance of phase contributions from the pulse envelope (5utt and the time- 
dependent phase — /3u(0() . Such a balance is provided by some value of pulse chirp. However, a sole phase 
balance is not sufficient as the pulse spreads. The spreading can be compensated by spectral filtering. Since the 
chirp causes the frequency deviation at pulse front and tale, the filter cuts off the higher- and lower-frequency 
wings of the pulse and, thereby, shortens it.^*''^^ 

The partial exact CSP solution of Eq. ([1]) is known for / (^l"*^!^) = k ^1 — <j (cubic-quintic nonlinear 

CGLE, for overview see^'^'^'*). In this work, the approximate method of integration of Eq. ([1]) in a general form 
will be proposed. The underlying approximations are 

Conjecture 2.1. T ^ y^, that is the adiabatic approximation; 

Conjecture 2.2. /3 3> a, that is the GDD prevails over the spectral dissipation. 

The first conjecture is valid for both solid-state and fiber oscillators operating in the PDR because the pulse 
is strongly stretched in the regime under consideration. The second conjecture is valid for both broadband solid- 
state (i.e. Ti:Sapphire^^ and Cr:YAG^^) and fiber oscillators.^^ Thin-disk oscillators based on the narrowband 
active media^^ can approach the limit of /3 « a and this issue will be addressed in Subsection 2.4. 

2.1 Cubic nonlinear CGLE 

The simplest version of ([T]) corresponds to x = and the SAM function is^^ 



Such a function approximates nonlinearity of low-energy solid-state and fiber oscillators. The K-parameter 
describes a nonlinear gain due to loss saturation. Eqs. (jll2p lead to a dissipative generalization of the nonlinear 
Schrodinger equation. 

Let's make the traveling wave reduction of Eqs. p|2p by means of ansatz 



2. DISSIPATIVE SOLITON OF NONLINEAR CGLE 




(1) 




(2) 



u (z, t) — \J P {t) exp [icj) {t) — iqz] , 



(3) 



where P is the z-independent instant power, (j) (t) is the time-dependent phase, and q is the phase due to shp 
of the carrier phase with respect to the envelope. In contrast to the tradition approach, we do not impose 
hereinafter any restriction on the time-dependence of phase. Substitution of ^ in ^ supplemented with ^ 
as well as taking into account the approximations under consideration and the obvious restrictions P > and 
< oo lead to 



jP{t) = /3A2(l-tanh2[<K(l-Fc)A/37]) 

aA^ = 7^(0) (4) 

2 — c 

where = q/P, P (0) is the CSP peak power, and the control parameter c = aj/f3K. One can see, that, on 
conditions that the appropriate normalizations are used, the CSP is two-parametric and depends on only a and 
c. 

Since a pulse is strongly chirped in the FDR, one may treat as a rapidly varying function and apply the 
method of stationary phase to the Fourier image of w.^^ As a result, the spectral power is 



piu;)^\e{.)f^j^^H{A^^u;'), (5) 



where e (w) = T[u] = J dty^ P (t) exp [i(f) (t) — iuit] and H (x) is the Heaviside function. That is the spectrum is 
flat-top and truncated at ±A. The latter parameter plays a role of the spectrum half-width. From Eq. ^ the 
CSF energy is 

Pdt= p=^^. (6) 
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Important features of ^ are i) < c < 2, ii) there is no physically nontrivial limit cr = 0, and iii) spectral 
width increases with a and c. Hereinafter, the solutions with the properties of ii) and iii) will be termed as the 
Schrodinger (or negative) branch of CSF. 

2.2 Cubic-quintic nonlinear CGLE (x = 0) 

Let's consider the case with negligible higher-order SFM and with the SAM in the form of^^ 



f[\u\')=^[l-<;\u\')\u\\ (7) 

Here the ^-parameter corresponds to the SAM saturation. Such a form of SAM can be attributed to the Kerr- lens 
mode locking or the mode-locking due to polarization modulator. The former technique uses power-dependence of 
the laser beam size due to self- focusing inside a nonlinear medium. As a result, the overlapping between the laser 
and pump beams becomes power-dependent, as well. This leads to the SAM in the form under consideration. For 
the mode-locking technique utilizing the power-dependent polarization, Eq. ([7]) can be considered as a low-order 
in P approximation of the trigonometric SAM function. The SFM saturation (i.e. the x-term) can be omitted if 
the beam confocal length inside an active crystal is much less than the crystal length (in a solid-state oscillator) 
or the mode is strongly confined (in a fiber oscillator). Also, x = for an airless thin-disk oscillator. 

Substitution of ([3]) in ([1]) supplemented with ([7]) and taking into account the approximation under consider- 
ation result in^^ 



7P = q~pn\ (8) 



where il = 0^ is the frequency deviation from the carrier frequency cuq. 
The restrictions P > and fit < oo allow obtaining 



7P(0) = /3A2 = ^(l-c/2±V(l-c/2)"-4a^/K 



a = ^(A2-172)(02 + S2), (9) 



= ^(l + c)-^7^(0). 

Integration of second Eq. ^ in the combination with first Eq. ^ gives the implicit expression for CSP profile.^® 
To define the spectral shape of CSP and its energy, one may use the approach described in previous Subsection. 
The spectral power is 

p (uj) ~ '-^^ 10 

That is the CSP spectrum has the truncated Lorentz profile, where A plays a role of the spectral half-width (if 
S > A, otherwise the spectral width is defined by S). 

Eq. PH)) allows obtaining the CSP energy 

= -^arctan . (11) 

Hence, Eqs. (IQIlOllip define the CSP completely. One can see, that the CSP parameters depend on only 
two control parameters: c = aj/^K and a = a^^/n. Since the saturated net-gain is energy-dependent due to 
gain saturation, one may choose c and E* as two control parameter, where E* is the energy of steady-state 
(continuous- wave, CW) solution of linearized Eq. ([T]). 

Thus, the CSP characteristics can be mapped on two-dimensional diagram (so-called "master diagram"), 
that makes theirs easily traceable (see below). 

2.3 Cubic-quintic nonlinear CGLE (x / 0) 

Eqs. (|1I7[) with x 7^ describe an oscillator operating under condition, that the beam is not confined strongly 
in a nonlinear element and its size variation affects the effective SPM. For instance, that can be a solid-state 
oscillator with the beam confocal length approaching or exceeding the crystal length; a high-energy oscillator 
with nonlinear plate providing the SPM; a fiber oscillator with no strong mode confinement, etc. 

The method sketched in previous Subsections leads to 

P = (12) 
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Figure 1. Positive branch: the CSP profiles Pit) (solid and dotted curves) and the corresponding chirps Q,t{t) (dashed 
curves and crosses) for the different values of quintic term b. a — 0.01, c — 1. 



where ^ (t) = y ^ + ^ ~ i^))' ^ = ^l/x ^^nd the following normalizations are used: t = tn^J K/a<^/i,, 

P = <,P, VL^ = ^a^jn, A ^ = A^a<r/K (the primes will be omitted below). Eqs. (fT2|) represent only solutions, 
which tend to those of Eqs. ^ when x — * 0. 

Expression for the dimensionless spectral power is 



TT (A - 1) {{A - 1) c6 + 4 {2u? - A2)) U (A2 - cj^) 
cA {[A - 1) (c (a + + 52 + cj2) + 6 (A2 - ^^2)) _ 2 (6 + 1) (A2 - t^2)) ' 



(13) 



where A = y 1 + ^^^75-^- One can see, that the CSP becomes three-parametric due to an appearance of 
non-zero x & 7^ 00). 

The CSPs under consideration subdivide into two classes accordingly two signs in Eqs. ()9I12|) : i) positive 
branch and ii) negative (or Schrodinger) branch. The negative branch can be transformed into the soliton- 
like solution of dissipative nonlinear Schrodinger equation, when the quintic nonlinear terms tend to zero (see 
Subsection 2.1 and^^^). The positive branch, possessing physically nontrivial limit a = 0, does not allow such a 
transformation. The additional characteristics of these two branches will be described in detail below. 

Let's consider the CSP profiles obtained from integration of (fT2|) . Figs. [T] and [2] show the profiles of the 
positive and negative branches, respectively. The 6-parameter is scalable and zero quintic SPM corresponds to 
h — !■ ±00. When 6 > 0, nonlinear phase shift increases with power. Such a sign of h corresponds, for instance, 
to a Kerr-lens mode-locked oscillator. One can see, that the peak power and, correspondingly, the chirp increases 
(decreases) with h for positive (negative) branch. When 6 < (saturable SPM), both flat-top (gray solid curve. 
Fig. [1]) and parabolic (dotted curve. Fig. [J) profiles appear. 

As it has been pointed (see Eqs. (|5I10I13|) ). the CSP spectra are truncated at some frequency ±A. The 
spectral profiles of positive branch are shown in Fig. [3l One can see, that there exist next spectral types: i) 
parabolic- top22' 27 (golid curve and circles; the last corresponds to large contribution of saturable SPM), ii) finger- 
like22'27 (dotted curve; that is the truncated Lorentz profile), and iii) concave2^' 2^' (dashed curve). In contrast 
to the model presented in Ref.,2^ our model predicts that the concave spectra, which are widely presented in the 
fiber oscillators, are stable. The last conclusion means that such spectra exist for the saturable SAM (i.e. for 
> 0). The cause of concave spectra is the positive quintic SPM (i.e. the SPM growing with P). 

Spectra of negative branch are shown in Fig. (4] The spectra are typically narrower than those of positive 
branch. Therefore, the finger-like profiles disappear and the typical shapes are i) parabolic- or flat-top, and ii) 
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Figure 2. Negative branch: the CSP profiles P{t) (soUd curves) and the corresponding chirps Q.t{t) (dashed curves) for 
the different values of quintic term b. a = 0.01, c = 1. 
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Figure 3. Positive branch: the CSP spectra for the different values of quintic term b. a = 0.01, c = 1. 
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Figure 4. Negative branch: te CSP spectra for the difFerent values of quintic term h. a = 0.01, c = 1. The profiles for 
b = 20 and -5 coincide. 
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Figure 5. The normalized spectral half-width in dependence on the saturated net-loss parameter a for the different values 
of quintic SPM. c = 1. Positive branch - solid curves; negative branch - dashed curves. 



concave. The latter shape correspond to wider spectra formed in the presence of positive quintic SPM (i.e. the 
SPM growing with P). 

Dependence of the normalized spectral half-width A on the normalized saturated net-loss parameter a for 
the different h is shown in Fig. [S] When 6>Oor6<— 2, the positive branch CSP has a wider spectrum than 
the negative branch one. Under these conditions, the spectrum width of positive (negative) branch decreases 
(increases) with the a-growth. The behavior of negative branch corresponds to that of CSP in the dissipative 
nonlinear Schrodinger equation (see Subsection 2.1). 

The CSP exists within confined region of a, which narrows with the decreasing positive h. When > 6 > —2, 
the positive branch exists only for nonzero a. This means that only negative-branch CSP can spontaneously 
develop in an oscillator with > 6 > — 2. 

Figs. I6|7I show the master diagrams, which are the two-dimensional slices of three-dimensional parametric 
space of CSP (for x 7^ 0). The dimensionless parameter E equals to the dimensional energy E multiplied by 




Figure 6. Master diagram: b — > oo (solid curves), b = 0.2 (dashed curves). Black curves are the stability thresholds 
(a = 0) , gray curves are the border between the positive (+) and negative ( — ) branches, light gray curves are the positive 
and negative branches for a = 0.01. 



(k/7) ^fn^Joi and can be easily related to E* by means of expansion of a in the vicinity of laser threshold:^^ 
a « 5{E/E* — 1) [5 = da /dE\^^^,). The black curves (the solid ones in both Figs, and the black dashed 
one in Fig. [6]) correspond to the stability threshold against the continuum excitation, that is a = along this 
curve. The stable CSP exists below these curves (i.e. the curves correspond to the maximum values of c for a 
given E). One can see, that the stability threshold c decreases with E. Physically, that means, for instance, 
the growth of GDD required for the CSP stabilization (c (x 1//9). It is important, that there is an asymptotic 
value of c (i.e. an extra-growth of E does not change the threshold value of c substantially). The contribution 
of positive (negative) quintic SPM narrows (broadens) the stability region. Also, the region, where the positive 
branch exists, narrows for the growing positive x with the decrease of positive b). It should be noted, that 
the positive branch disappears for small a, when > 5 > —2. 

The master diagrams demonstrate main difference between the CSP branches. When 6 > is not too small, 
the asymptotic behavior of isogains (i.e. the curves of constant a) with the i?-growth demonstrates that the CSP 
energy is scalable. That is E cx E* and the proportionality coefficient is weakly dependent on the parameters 
of ([1]). In this sense, the negative (Schrodinger) branch is not energy-scalable, because the energy depends on c 
weakly and the change of E requires a substantial change of c (e.g. a substantial GDD growth as c (x 1//3). 

2.4 Generalized nonlinear CGLE 

In this subsection the SAM corresponding to a perfectly saturable absorber^^' '^'^ will be considered. This type 
of SAM represents a semiconductor saturable absorber mirror (SESAM), which are extensively used in CPO 
oscillators. ^^"^^ If the CSP width (usually, few picoseconds) excesses the SESAM relaxation time T,. (hundreds 
of femtoseconds), the SAM fimction can be written in the form 



+ ^\u\ 



2 ' 



(14) 



where k = /i<;, fj. is the modulation depth, — Tr/EgS is the inverse saturation power [E^ is the SESAM 
saturation energy fluency, S is the beam area on SESAM). We assume below that x 

Eqs. (|lll4p can be reduced by the above described method to 
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Figure 7. Master diagram: b = —5 (solid curves), b — —1.8 (labeled gray curve). Black curve is the stability threshold 
(a = 0) , gray curve is the border between the positive (+) and negative ( — ) branches, light gray curves are the positive 
and negative branches for a = 0.01. Gray curve b = —1.8 corresponds to the negative branch. 
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Here a = <t/ii, c = aj/ Pk, T = (2 — c) 
and the CSP energy are 



— 4a (2 — a + c). The spectral profile (truncated parabolic- or flat-top) 
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Using the normalizations from previous Subsection reduces Eq. (|17p to 
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where primes for the normalized values are omitted. 

The master diagram following from Eq. (fT8|) is presented in Fig. [H The structure of diagram is similar to 
that described in previous Subsection, but there are two important differences: i) c-parameter asymptotically 
tends to zero with growth of E, ii) negative branch is not energy-scalable within all range of existence. The 
latter conclusion is obvious from behavior of isogains in Fig. [5] as one can see that the c-scaling weakly affects 
E. Thus, the energy remains almost constant along this isogain when the GDD scales (c oc 1//3). However, the 
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Figure 8. Master diagram of generalized nonlinear CGLE. Black curve is the stability threshold (a = 0) , gray curve is 
the border between the positive (+) and negative (— ) branches, light gray curves are the positive and negative branches 
for a = 0.01 and 0.1. 
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Figure 9. Normalized spectral half-width in dependence on c for the different values of a. Positive branch - solid curves; 
negative branch - dashed curves. 

energy scaling for the negative branch can be provided by a simultaneous growth of a and (3 so that c remains 
constant (see also the normalization for E). 

For the negative branch, the spectrum narrows with the c-decrease due to growth of the GDD contribution 
(Fig. [9|), which stretches the pulse when the energy remains almost constant. When E changes weakly along the 
isogain corresponding to the negative branch, the spectrum broadens with a (Fig. [51 c oc a). The explanation 
is that the growth of spectral filtering enhances the cutoff of red (blue)-shifted spectral components located on 
the pulse front (tail). The growth of cutoff shortens the CSP and, for a fixed energy, P (0) increases. Since 

oc P(0), the spectrum broadens. For the positive branch, the spectrum initially broadens with c-decrease 
(Fig. [9]). This can be explained as a result of i?-increase, which is necessary for keeping in the isogain (Fig. [8]). 
That is, the SPM contribution increases and the spectrum broadens. However, further decrease of c narrows the 
spectrum due to growth of the CSP width. The latter results from either GDD growth (c (x 1//3) or suppression 
of spectral cutoff on the pulse wings due to a-decrease (c cx a). 



The CSP spectral width behaves similarly for the cubic-quintic nonlinear CGLE with x = 0. It should be 
noted additionally, that the positive branch remains energy-scalable for the SAM under consideration, but there 
is a need in the decrease of c to provide such a scaling (Fig. [5]). 

One has to emphasize that the master diagrams cover both solid-state and fiber oscillators. In a fiber 
oscillator, the GDD value is larger in comparison with that in a solid-state CPO, but the spectral filter bandwith 
(~25 nm) is small. As a result, the excess of the ratio /3/a over that for a Tiisapphire CPO is only tenfold.'^" 
Simultaneously, an excess of the ratio 7/K over that for a Tirsapphire CPO is tenfold as well. As a result, the 
c-parameter is ~ 1 and, dynamically, there is no substantial distinction in kind between the ANDi fiber and 
the solid-state CPOs. One difference is that a larger GDD and a comparatively smaller E shift the operational 
point of a fiber oscillator into the negative branch region. While the operational point of a solid-state oscillator 
belongs to the positive branch. As a result, their scaling properties are different. 

In the case of narrow-band thin disk solid-state oscillator, the assumption /3 3> a is violated. The result is 
the smoothed spectrum edges. Such a smoothing increases with the decrease of SPM (e.g., when a resonator 
becomes airless). Nevertheless, the numerical simulations demonstrate that the analytical model describes the 
spectral width of CSP adequately even for a ~ /?. However, the case of a > /3 is beyond the bounds of the theory 
under consideration. 

3. CONCLUSIONS 

The completely analytical theory of CPOs has been developed. It has been found, that such oscillators (both 
solid-state and fiber) can be modeled by the nonlinear CGLE, which has been integrated approximately on the 
basis of the proposed method. As a result, one may easily trace the properties of CSPs, because the CSP has 
two- (maximum three-) dimensional representation (so-called, master diagram). Four types of SAM have been 
analyzed on the common basis. These SAMs are inherent in both solid-state and fiber CPOs mode-locked by 
either self- focusing, or SESAM, or polarization modulator. The CSPs have been found to be subdivided into 
two classes with different scaling properties. This properties have been traced and described within the whole 
parametric space. 



The mathematical apparatus is presented in detail in 


http : / /info . tuwien . ac . at /kalashnikov /NCGLEl.html 


and 


http: / /info. tuwien.ac.at/kalashnikov/genNCGLE. html 
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